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Abstract 

We present a supersymmetric extension of the Standard Model (USSM-A) witli an anomalous 
U{1) and Stiickelberg axions for anomaly cancellation, generalizing similar non-supersymmetric 
constructions. The model, built by a bottom-up approach, is expected to capture the low-energy 
supersymmetric description of axionic symmetries in theories with gauged anomalous abelian inter- 
actions, previously explored in the non-supersymmetric case for scenarios with intersecting branes. 
The choice of a USSM-like superpotential, with one extra singlet superfield and an extra abelian 
symmetry, allows a physical axion-like particle in the spectrum. We describe some general features 
of this construction and in particular the modification of the dark-matter sector which involves 
both the axion and several neutralinos with an axino component. The axion is expected to be 
very light in the absence of phases in the superpotential but could acquire a mass which can also 
be in the few GeV range or larger. In particular, the gauging of the anomalous symmetry allows 
independent mass/coupling interaction to the gauge fields of this particle, a feature which is absent 
in traditional (invisible) axion models. We comment on the general implications of our study for 
the signature of moduli from string theory due to the presence of these anomalous symmetries. 
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1 Introduction 



Extensions of the Standard Model (SM) describing axion-like particles - and with supersymmetry as a 
basic low energy scenario - are an interesting area of investigation which has the potentiality to provide 
an answer to a series of unsolved theoretical issues. Among them are those concerning the possible 
presence of anomalous extra neutral gauge interactions at current and future colliders in some special 
channels, especially in the search for an anomalous extra Z' . This investigation could also clarify the 
role of weakly coupled pseudoscalars in the early universe. For this reason several studies addressing 
the experimental detection of pseudoscalars at future experiments [D El O [H [5l [6] , has received an 
impressive impulse in the recent literature. 

One of the distinctive features of these extensions is the presence of extra abelian interactions 
which are anomalous. We just recall that anomalous C/(l)'s are quite common in several string 
constructions and that the mechanism of anomaly cancellation, if realized at low energy by a Wess- 
Zumino counterterm (WZ), may cause the presence of a physical axion in the spectrum. This result 
points directly towards the possibility of having a new dark matter candidate (see also [7j), which is 
certainly one of the most appealing features of this class of theories [8] . 

One of the first successful realization of the non-supersymmetric version of these models comes 
from special vacua of string/brane theory (orientifold vacua), in the form of stacks of intersecting 
branes, which induce a gauge structure given by the product of U{N) ~ SU{N) x C/(l) factors, where 
N is the number of branes of each stack (see [9] for an overview). Among the U{1) factors, one of 
them is identified with the SM hypercharge (U{1)y), while the remaining ones are anomalous and 
involve Stiickelberg axions for anomaly cancellation. In effective string models the abelian structure 
is in general characterized by the presence of several U{1) factors, described in the hypercharge basis 
by direct products of the form Gi = U{1)y x U{l)i x ... x U{l)p, with an anomaly-free hypercharge 
generator and p anomalous ?7(l)'s which are accompanied by axions bi, with i = 1,2,... p. The 
anomalous ?7(l)'s in this construction are in a broken phase, called the "Stiickelberg phase". In 
particular, after electroweak symmetry breaking (EWSB), one of the axions becomes physical ^ and 
is characterized by independent mass/coupling relations, where the coupling appears in an ordinary 
bFF interaction with the gauge fields, providing a generalization of the Peccei-Quinn (PQ) axion. One 
shortcoming of this description, at this time, is the absence of a supersymmetric extension of it with 
the appearance of a physical axion. The generalization to the supersymmetric case of these theories is 
interesting on several grounds. For instance, it allows to study an entire new class of extensions of the 
MSSM in the presence of a gauging of the axionic symmetries [10] and, at the same time, represents an 
intermediate step toward the unification with gravity of the same models, within certain formulations 
of supergravity |1H I12j. The formulation of [TU], which is specific for a MSSM superpotential parallels 
a previous general study of the same topic contained in [TT| . 

Therefore, these types of constructions provide a consistent framework for the study of the effects 
of moduli of string/brane theory within scenarios with large extra dimensions or via supergravity. 
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together with their low energy imphcations in cosmology and in collider physics [13j. Recently, an 
extension of the MSSM containing an anomalous U{1), made massive by a Stiickelberg supermultiplet 
|14j has been introduced in [lOj. This has been based on the superpotential of the MSSM with an extra 
abelian symmetry. One of the features of this construction is the absence of a Higgs-axion mixing, since 
the bosonic component of the Stiickelberg multiplet remains an ordinary goldstone mode. Therefore, 
the final theory is characterized by a physical axino but not by a physical axion. The objective of 
our analysis is to show that a similar construction can be performed in more general ways, thereby 
generating a model with a physical axion-like particle. This provides a complete supersymmetric 
generalization of the (gauged) PQ axion. We will work out the requirements that are needed in order 
to make this extension possible, detailing some of the arguments that have been presented in short 
form in [15] and analyzing the main features of the effective action of such a theory, that we call the 
USSM-A due to the anomalous U{1) (A) and to the specific choice of the USSM superpotential. 

Our work is organized as follows. We briefly describe the class of models that we are going to 
investigate, outlining their basic structure, together with their supersymmetric generalizations. Along 
the way, we will underline the differences between our construction and the previous construction of 
[lOj . We show how a physical axion is bound to appear in the spectrum and describe all the sectors 
of this theory. We derive the corresponding generalized Ward identities and characterize the Chern- 
Simons interactions of this class of models bringing up one typical example of application. We study 
the neutralino sector of the model and present a brief numerical analysis of its spectrum. Most of our 
attention in this work focuses on the basic characterization of this model, stressing on the mechanism 
that allows a physical axion in the spectrum. We conclude with some comments on possible extensions 
of this analysis to more general potentials characterized by moduli in different scenarios derived from 
string theory. 

2 Supersymmetric Extensions of the Standard Model and extra 

[/(l)'s 

Abelian (anomaly-free) supersymmetric extensions of the SM have been discussed in several previous 
works [16t ll7 t [TS lll9t [ ^ l21|. In [2D] the authors explore an extension of the Minimal Supersymmetric 
SM (MSSM) with an extra SM singlet chiral superfield S, with chiral charges chosen so to allow 
trilinear couplings of S to the two Higgs doublets Hi,H2 in the superpotential. The /i term, in this 
case, is generated by the vev of the scalar component of 5, precisely by the SHi ■ H2 interaction. 
The structure of this model, usually called USSM, shares some similarity with the nearly-Minimal 
Supersymmetric SM (nMSSM) [22] and the next-to-Minimal Supersymmetric SM (NMSSM) [23]. In 
all of these three models the extra scalar S is introduced for the same purpose but in the nMSSM and 
NMSSM this field is a singlet under the complete gauge group (which is the same as the SM) while 
in the USSM the field is charged under the extra C/(l). We recall that the nMSSM and the NMSSM 
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differ at tlie level of the superpotential in the structure of the pure S contribution, which is either 
linear (nMSSM) or cubic (NMSSM). 

In the approach of [20] this appears to be a necessary requirement since a scalar superfield, singlet 
under the complete SU (3) x SU (2) xU{1)y xU{1)b gauge group, while solving the ^ problem, however, 
does not allow a consistent pattern of EWSB, leaving the extra Z' of the neutral sector massless. This 
construction is realized with an anomaly-free chiral spectrum. 

2.1 MSSM and USSM with an anomalous f/(l)' 

In [To] the authors investigate a supersymmetric extension of the SM with an extra U{1), based on 
the superpotential of the MSSM. They make an important step forward in the analysis of this class of 
theories, using a bottom-up approach, that is by 1) fixing the effective action of the anomalous abelian 
symmetry using the Stiickelberg supermultiplet to give mass to the anomalous gauge boson and 2 ) 
using Wess-Zumino counterterms to balance the mixed and cubic ^7(1)^ anomalies of the theory. A 
third element of the construction is the possible presence of Chern-Simons interactions [8] which find 
their way to low energy from string theory [24j , and which amount to a re-distribution of the anomaly 
starting from a symmetric distribution on each leg of the anomaly vertex. This re-distribution is 
allowed whenever the symmetry of the vertex does not allow to uniquely define the breaking of the 
Ward identities separately on each of its legs. The meaning of this freedom, from the point of view of 
the effective field theory, is that each model allows a set of additional (defining) Ward identities for 
the distribution of the anomaly which are a specific feature of anomalous models in which the trilinear 
gauge interactions are not identically zero (in the massless fermion phase, the chiral phase). 

In the first supersymmetric version of these models |10j . the ordinary MSSM Lagrangean is nat- 
urally extended by the Stiickelberg multiplet which provides a kinetic term for the same multiplet 
while rendering the extra Z' massive. The defining phase of the model is, therefore, the Stiickelberg 
phase. In this construction the bosonic partner of the axino, which is the fermionic component of the 
multiplet, remains a goldstone mode after EWSB and is therefore unphysical. 

2.2 Inducing Higgs-axion mixing 

At the origin of the physical axion is the mechanism of Higgs-axion mixing. For this to take place one 
needs a Higgs sector which is charged under the anomalous U{1)b so that the mass of the anomalous 
gauge boson comes from a combination of the Higgs and Stiickelberg mechanisms. In the case of the 
MSSM this mixing does not occur even if the two Higgses are charged under the anomalous U(l). 
The presence of a /i term in the superpotential forces the two charges of the two Higgs doublets to 
take opposite values, thereby guaranteeing also the cancellation of the extra anomalies due to the 
circulating higgsinos, but is not enough to give mass to the anomalous gauge boson. In other words, 
in the absence of a Stiickelberg multiplet the mass matrix of the gauge boson has still an additional 
null eigenvalue. The true mechanism of mass generation of the anomalous Z', therefore, is just the 
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Stiickelberg, which in this situation is a goldstone mode. In fact, one reobtains a massive Yang- Mills 
theory just by going to the unitary gauge and eliminating the axion. 

3 The structure of the model 

A simple way out in order to have Higgs-axion mixing and a light axion in the physical spectrum 
consists to use a modified superpotential as in [20] , but now with an anomalous gauge structure, and 
to combine it with the Lagrangean of the Stiickelberg supermultiplet. In other words, we move from 
the superpotential of MSSM-type to the one typical of the USSM, introducing an extra scalar superfield 
S which is non-singlet under an extra U{1)b, maintaining the anomalous structure induced by the 
extra neutral current. This specific assumption allows to remove the second massless eigenvalue in the 
mass matrix of the gauge bosons and allows to induce Higgs-axion mixing once that the Stiickelberg 
mechanism is invoked to contribute to the mass of the extra Z' . The conditions that need to be 
verified in order to have a physical axion in the spectrum are obtained from an analysis of the CP- 
odd sector of the theory and involve both the potential and the derivative couplings (mixings) of the 
massive gauge bosons with their goldstones [ZidGzi) extracted from the broken phase. In general, the 
presence of extra singlet superfields in the superpotential allows such a mixing and we will illustrate 
this requirement in one of the sections below. The analysis that we will present in the next sections 
has the goal to clarify this point, starting from the MSSM case, where none of the CP-odd states 
acquires an axion-like coupling. 

These new features do not affect the chargino sector with respect to the MSSM. 

4 The superpotential 

The construction of models characterized by a physical axion in their spectrum requires an appropriate 
superpotential. In order to obtain this, we consider the introduction of an extra SM singlet S. For 
this reason, the superpotential of the model investigated is given by 

W = XSHi-H2 + yeHi-LR + ydHi-QDR + yuH2-QUR, (1) 

which coincides with the model of ^0|, called the USSM. We refer to Tabled] for a list of the charge 
assignment of the chiral superfields of our model; the scalar superfields corresponding to SU (3), SU{2), 
U{1)y and U{1)b are, respectively, G"^(x,6',(9) (with a=l,2. ..,8), W'{x,e,9) (with i=l,2,3),y (x, 6*, ^) 
and B{x, 6, 6) and they fall in the usual adjoint representations of the gauge group factors. 

We have denoted the charges by Qf,x, where X denotes the hypercharge (Y), the charged 
bosons (W), the non abelian gluons (G) and the anomalous gauge boson (B). At the same time we 
denote with Bx the charges of the X superfield respect to the anomalous U{1). Unlike the NMSSM 
and the nMSSM, W does not contain linear and cubic terms in S in order to preserve the gauge 
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Superfields 


SU(3) 


SU(2) 


U{1)y 


U{1)b 


h{x,e,e) 


1 


1 





— 


s{x,9,e) 


1 


1 





Bs 


L{x,e,e) 


1 


2 


-1/2 


Bl 


R{x,9,e) 


1 


1 


1 


Br 


Q{x,e,e) 


3 


2 


1/6 


Bq 


Unix, 9,6) 


3 


1 


-2/3 


Bur 


DR{x,d,e) 


3 


1 


+1/3 


Bdr 


Hi{x,e,d) 


1 


2 


-1/2 


Bhi 


H2{x,e,e) 


1 


2 


1/2 


Bh2 



Table 1: Charge assignment of the model; boldface numbers indicate the dimensions of the corre- 
sponding representations. 

invariance in the presence of a non vanishing Bs charge. This requirement is strictly necessary if the 
extra scalar S is only a SM singlet. Gauge invariance gives the conditions 

Bh,+Bh,+Bs = 

Bh,+Bl + Br = 

Bh,+Bq + Bdr = 

Bh,+Bq + Bu^^ = 0, (2) 

which will be used below. It is not hard to show that the possibility of declaring S" to be a singlet 
under the entire gauge group {Bs = 0) SU{3) x SU{2) x Gi leaves an extra gauge boson massless 
beside the photon, after EWSB and as such it is not acceptable. 

4.1 Anomaly cancellation: defining the model 

We start by identifying the anomalous contributions of the model, whose gauge structure is of the 
form SU{3) x SU{2) x C/(l)y x U{1)b. 

The anomalous trilinear gauge interactions are all the ones involving the extra anomalous U{1)b, 
namely {U{1)b, U{1)b, U{1)b}, {U{1)b, U{1)y, U{1)y}, {U{1)b, U{1)b, U{1)y}, {U{1)b, SU{2), SU{2)}, 
{U{1)b, SU{3), SU{3)}. In terms of the charges we can write each sector as follows 

■Abbb = ^ Q),b 
f 

Abyy = X] Qf,B Q},y 

f 

Abby = ^ Q\b Qf,Y 

f 
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/ 

^BGG = 5^Q/,BTr [t"T^] , (3) 
/ 

where are the generators of SU{3) and r* the PauU matrices. Compared to the analysis of |10j . 
here we have anomalous trilinear interactions also in the sector involving the SU{3) mixed anomaly 
due to the non vanishing charge Bs- Using the constraints coming from the Yukawa couplings and 
the conditions of gauge invariance, the expressions of the anomalies take the form 

Abbb = ?>{QBl+Wfj^ + Wl^) + {QBl+3Bl) + {2Bl^+2B^H,+Bl) 

= -3B^H, - K'^Bl + 18Bq - 7Bs)Bl^ - 3{3BI + (185q - 7Bs)Bs)Bh, 
+3Bl + Bs{27Bl - 27 BsBq + 8B|) 
Abyy = 3{6BqY^ + 3Bu^yS^ + 3Bd^Y^J + {6BlYI + 3BrY^) 

+{2Bh,+2Bh,)Y^, 

= ^{-3Bl-9Bq + 7Bs) 
Abby = KQBIYq + 3Bl^Yu^ + 3BI^Yd^) + [QBIYl + 3BIYr) 

+{2Bh^-2Bl^)YH, 

= 2BhA3Bl + 9Bq - 5Bs) + {12Bq - 5Bs)Bs 
Abww = \{ISBq + QBl + 2Bh,+2Bh,) = 3Bl + 9Bq-Bs 

Abgg = \{QBq+3Bu^ + 3Bd^) = ^Bs, (4) 

where Yq, Yl are the hypercharges of the left-handed doublets of the quarks and leptons respectively, 
while Yu^^ Yd^, are the hypercharges of the fJ/j, Dr^ R superfields which correspond to the hyper- 
charges of the right-handed quarks and leptons, with the opposite sign. 

In the absence of a specific charge assignment coming from a string (or other) construction, these 
equations can be interpreted as defining conditions of a specific model. The role of string theory or of 
any other construction is to fix the charges, but for the rest the basic structure remains determined 
by the approach outlined below, and as such is truly general. 



5 The Stiickelberg multiplet 

In supersymmetric models the cancellation of the anomaly using the Wess-Zumino (WZ) counterterm 
can be obtained by the introduction of a Stiickelberg supermultiplet, associated with the extra U{1). 
The multiplet contributes to the supersymmetric version of the Stiickelberg mechanism and in 
the WZ interaction that describes the coupling of the supermultiplet to the gauge supermultiplet. We 
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recall that in anomaly-free theories the Stiickelberg mechanism has the feature of contributing to the 
mass of the anomalous gauge boson, eventually also in combination with the Higgs sector |25l 126] . 
This construction holds both in the non-supersymmetric and in the supersymmetric case. 

Obviously, the presence of a mixing between the Higgs and Stiickelberg components in the potential 
of more generic models in an anomaly-free theory, produces a new CP-odd component in the scalar 
sector, but deprived of axion-like couplings. On the contrary, these couplings appear in the case in 
which the two mechanisms (the Higgs and the Stiickelberg) involve an anomalous U{1), due to the 
presence of Wess-Zumino terms, for specific superpotentials. These interactions are induced in the 
effective action by the mechanism of anomaly cancellation. 

The Lagrangean describing the Stiickelberg supermultiplet is given by 



£., = /^.[2M.,B + b + (,f (5) 

where B is the abelian scalar superfield associated to the extra C/(1)b, b is a left-chiral superfield and 
Mst is the Stiickelberg mass. 

The former Lagrangean is invariant under the following gauge transformations 

B B' + i(^A-A^' 

h h'-i2Mstk (6) 
where A is a generic left-chiral superfield. Introducing the component fields expansion we obtain 

B = -ea''eBf, + ieed\B-imXB + ^9999DB (7) 



b = b + iV29iJb-i9(j''9d^b + —999a''df,'il^k--9999n\b + 99FY„ (8) 



Pi - 1 
—999a^dai)\, - 

where i?^ is the Stiickelberg field, Xb,^b are respectively the left- and right-handed Stiickelberg 
gauginos, Db is the corresponding D-teim for the gauge supermultiplet of B^ , 6 is a complex scalar 
field, V'b is the supersymmetric axion (axino) and Fb is the F-term of b. 

After the integration over the Grassman variables the Lagrangean density is given by 

Cst = 2{d^lmb + MstB^f + iVb^'^S^V^b + iV^b^'^S^V'b + 2F^fI + AM^t Re 6 Db 

-2V2Mst (VbAfi + h.c.) , (9) 

where the auxiliary fields Fb and Db will be defined in the next sections. 
5.1 The axion Lagrangean 

The axion Lagrangean contains the Stiickelberg gauge-invariant terms introduced above and the Wess- 
Zumino interactions for the anomaly cancellation and it is given by 



Caxion = \ j d''9{h + bt + 2M,tBf -\l d''9 



h)G Tr(gg)b + h)w Tx{WW)h 
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+6ybWjTy^'" + bnhW^W^''' + byB^W^W^^'' 6{9^) + /i.e.} , 

(10) 

where we have denoted with Q the field-strength of SU{3)c, with W the super symmetric field-strength 
of SU{2), with and with the supersymmetric field-strength of U{1)y and U{1)b respectively. 
The factors in front of the WZ counterterms (bx) are determined by the standard conditions of 
anomaly cancellation. The Lagrangean, in our case, contains extra WZ counterterms respect to |10j . 
in particular we need to impose the cancellation of the mixed B — SU (3) — SU{3) anomaly, which 
is now non-vanishing due to the charges of the two higgsinos in the model, which are not opposite. 
In the MSSM this cancellation is identical, due to the specific color charges of the fermions in each 
generation. This implies that in our case the effective action contains both a bGG interaction of the 
axion with the gluons and a vertex involving the corresponding gauginos (gluinos). 
Expanding the Caxion in the component fields we obtain 

Caxion = \ [d^lm b + MstB^f + ^VbCT^^^^b + ^V^^^^^/xV-b + - ^(V'bAi? + h.c.) 

-h^Ge^'f-^Gl^G-^^i-nb- Ibw e'^'^^'^Wi^W^Jmb - ^bye^^'"'"^ F^^Fjjmb 
-^tBe^-P^F^^F^^lmb- \byBe'^''^- F^^F^^Jmb 



+^6G[Jlm6 Xgao^D^Xga - ^^b Xg-o^o'^Gl^ + Ag. X^a + ^V'b A^aZ)" + h.c\ 

+6y[Im6Aya'^Z)^Ay - -l^V'bAya'^a'^Fj, + ^FbAyAy + -L^bAyDy + h.c\ 

- i 11 

+bB[lmbXB(7''Df,XB - ^VbABCT'^a^'F^^^ + -FbA^As + -^^j^XbOb + h.c] 

+bYB[(imbXYa''df,XB + ^Fy^XyXb + -^VbAyi^s - ^Aya'^^'^F^^^Vb) 

+ {¥ ^ B) + h.c], (11) 

with the F and D terms given by 



Fh = —{bcXgaXga + 6vkAvk»Avk* + fey Ay Ay -I- bBXBXB + ^y^AyAs), 
Db = -[^{BlL^L + BrR^R + BqQ^Q + BuUIUr + BdD^^Dr 

+Bh,h\Hi + BhMH2 + BsS^S) + ^MbsXB + &yi?Ay)], 

Dy = -[^(^^^ - 2^^^ - Iq^Q + P^Ur - Id^Dr + hIh, - HIH2) 

+^V'b(byAy + bYBXB)] 
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= -\[gs{Q^T''Q + uIt'^Ur + d^j.t'^Dr) + ^i^^Xga], (12) 

in which we have terms coming both from Caxion and from the USSM Lagrangean that can be found 
in the appendix. 

5.2 The kinetic mixing 

In these type of supersymmetric models the extra U{1)b sector can mix with C/(l)y in different ways. 
In particular, in the context of USSM — A, the kinetic mixing is treated as in the NMSSM with the 
inclusion of an anomalous C/(1)b symmetry and the extra singlet S is charged under B. 

The lagrangean for the gauge fields is modified by introducing a mixing term B — Y proportional 
to a small parameter sin a 

Cmixing = -\ j (fO 2sma W^'^W!^ 6\e) + h.c. (13) 

where sin a represents the mixing between the two abelian structures U{1)y and U{1)b- In the same 
way, the gauge mass terms lagrangean in the presence of kinetic mixing is modified by the inclusion 
of a term proportional to the mass parameter Myb as follows 

jC-GMTmix = \jd^^ [MybW^^W^ + h.c] 5\e , 9). (14) 

Furthermore, the USSM — A is affected by another source of kinetic mixing coming from the mixed 
counterterm proportional to byB in the expression of Caxion- Expanding this expression in component 
fields we observe that the multiplet b contains the complex scalar field b whose real part can be 
Re 6 7^ and it generates a kinetic mixing proportional to oc byB^^b gy gB, where the coefficient byB 
fixed by the anomaly cancellation procedure, goes like the inverse of the Stiickelberg mass and can be 
neglected in this first analysis (see Ref. [TO]). In our formulation we assume sin a = for simplicity 
and we will give a more detailed analysis of the kinetic mixing in the context of the USSM-A in a 
forthcoming paper [27]. 

5.3 The Fayet-Iliopoulos terms 

To be as more general as possible, in theories with U{l)s gauge superfields we should add to the 
lagrangean the following Fayet-Iliopoulos (FI) term 

CFI = ^yDy+^BDB. (15) 

which is allowed by symmetry reasons. Here (,y,^B are two coefficients, while Dy and Db are the 
D-terms corresponding to the U{l)y and U{1)b symmetry respectively. In our analysis we omit these 
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contributions even if a quadratically divergent FI always appears in a field theory at one loop [28]. 
The reason resides in the fact that, in the low-energy lagrangean there should be a counterterm, which 
compensates precisely both the divergent and the finite part of the one loop contributions (see Ref. 
|10j). We are also omitting the terms responsible for the cancellation of gravitational anomalies. A 
more comprehensive description will be given in |27j . 



Some of the notations used in our analysis are recalled in the appendix, here we just mention that 
the scalars of the model are denoted, as usual, by a tilde ( ~ ). It is convenient to combine the axion 
sector and the F and D terms extracted from the other sectors of the total Lagrangean of the model. 
This combination is in general defined to be the auxiliary Lagrangean, or Caux, which is given by 



-y, 



-\XS\\hIH2 + h\Hi) - vId^^DrQ^Q - ylVLRR^ 
-vIUrUIQ^Q - \yu [sQ^Hitjl + /i.e.) - \yd {SQ^H2D^^ + h.c:^ 

Aye {sVH2R^ + h.c}j - VdVu (uIhIHiDr + /i.e.) - yeVd (oiQ^LR + /i.e.) 

-\[9s{Q^T-Q + U^T'^Ur + D^j.T'^Dr) + ^iJbXgaf 



(LtZ - 2R^R - -Q^Q + -U]^Ur - -D^^Dr + hIHi - HIH2) 



'2^/2 2,^ 3 " 3 

+^^b(6yAy + hYBXB)? - [^{BlL^L + BrR^R + BqQ^Q + BuUjiUR 

+BdD^rDr + BhMHi + BhMH2 + BsS^S) + ^MbnXB + ^ysAy)]' 

+ ^[V'bV'b(&G W + bw^W'^W + (by + feyB)AyAy 

+b%\B\B + (by + bB)byBXyXB + bBbyB^B^y) + h.c] (16) 
where the expressions of the D terms are now determined by Eq. (jl2p . 



6 Goldstones of the potential and of the massive gauge bosons 

The identification of the goldstone modes of the model requires a combined analysis of the potential 
and of the bilinear mixing terms Z^dGzi ^ the broken (massive) gauge bosons. Naturally, the 
expansion near the vacuum is consistent if the stability conditions of the potential near the expansion 
point are satisfied. The neutral goldstone modes corresponding to the physical neutral gauge bosons 
after the breaking are part of the CP-odd sector together with other physical components, spanning 
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together the entire CP-odd space. In general, in this sector, the potential contains a set of "fiat 
directions" , which appear as goldstone modes of the matrix of its second derivatives. These goldstone 
modes do not necessarily coincide with the goldstone modes {Gz') identified from the bilinear mixings. 
This turns out to be the case if the Stiickelberg decouples from the scalar potential while it gives mass 
to one of the anomalous gauge bosons. To clarify this point it is convenient to move back to the 
non-supersymmetric case. 

The allowed structure of the potential involves 6-indcpcndcnt (V) and b- dependent (V) terms, 
just on the basis of the symmetries of the Lagrangean, given by 

V=Y, {f^lHlHa + XaaiHlHaf) - 2Ai2(i^Ii/l)(i^|i/2) + 2>!^^\H^ T2Hi\\ (17) 
a=l,2 

and 

+A2 [hIH2^ H\Hie''^'^'^'^"^'^^^ + A3 {h\h^ HlHie~'^'^'^^~'^^^^ + c.c. (18) 

respectively, where the sum over / is a sum over the Stiickelberg axions of the (several) anomalous 
[/(l)'s. In the supersymmetric case this second contribution is, in general, not allowed, although it 
might appear after supcrsymmetry breaking. This second term or "phase-dependent term" is directly 
responsible for Higgs-axion mixing and for producing a massive axion. The interesting point is that 
in the supersymmetric case (with 6 a real field), even if V is not allowed, we may still, under some 
particular conditions, end up with a physical axion in the spectrum, as we are now going to elaborate. 

As we have mentioned, the identification of the goldstones of the theory is necessarily done using 
the kinetic term of the scalars, including the Stiickelberg, which in this case takes the form 

Pm^iI' + \'Df,H2\^ + ^{d^h + MstB^f. (19) 

The expansion of this equation near the stable vacuum gives the usual bilinear mixings characterizing 
the derivative couplings of the physical massive gauge bosons to the corresponding goldstones; rather 
straightforwardly one obtains the combination 

MzZ^'d^Gz + Mz'Z"'d^Gz' + ... (20) 

with Gz and Gz' being the true goldstone modes of the theory. Notice, if not obvious, that while 
Gz is just expressed as a linear combination of the two CP-odd components of the Higgs, Gz' on the 
other hand takes a contribution also from 6, due to the Stiickelberg mass term. Therefore, one of the 
special features of the combination of the Higgs and Stiickelberg mechanisms is that in some cases the 
potential of the model - F is an example of this situation, since it does not not include a b field - is 
not sufficient to identify all the goldstone modes. Clearly, if both V and V are present, then Gz and 
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Gz' can be extracted from the total potential and coincide with the goldstone modes extracted from 
the bilinear mixings of (|19p and (|2Up . In this case the physical axion tmms out to be massive. We 
recall that the quadratic part of the CP-odd potential takes the general form 



VcP-odd = iluiHlluiHlb) M 



Imi/0 I (21) 



for a suitable M matrix whose explicit expression is important but not necessary for our dicussion. In 
the case of the MSSM the structure of the potential coincides with that of V and one identifies only 
one physical CP-odd Higgs (called in the MSSM) which will not have an axion-like coupling, as 
can be verified by also a simple counting of the degrees of freedom before and after EWSB. In this 
case the orthogonal transformation that diagonalizes the CP-odd scalar sector takes the form 

' ' ^ O 22 



and involves the physical (massive) CP-odd Higgs A^ and a golstone mode G^. The above discussion 
goes through in a similar way also for the anomalous U{1) extension of the MSSM discussed in |lUj . 
For the case of a potential such as Vcp-odd = V + V' instead, there is indeed a mixing between the 
components of the Higgs and b and the diagonalization of the quadratic part of the potential gives 

j^Imi/0 j = 03 I^G? j (23) 

with O3 being an orthogonal matrix. We have denoted the physical field by x the NG-bosons 
by G?2- In this case it is rather obvious that x acquires an axion-like coupling, inherited from b. In 
other words h has an expansion in terms of x, G^ and G2 or, equivalently, in terms of Gz and Gz', 
where Gz and Gz' are identified by Eq. ([2D|) . The decomposition is clearly gauge dependent. One 
important comment concerns the nature of the bFF interactions in this case. 

In the unitary gauge the only axion-like couplings left involve the physical component of b, denoted 
by X) called "the axi-Higgs", which gives typical x^^F interactions. As we have mentioned above, in 
the absence of y , b decouples from the rest of the Higgs sector in M. In this case in the unitary gauge 
all the anomalous couplings can be removed, and the theory goes back again to its original anomalous 
form, with the old Lagrangean now replaced by an ordinary massive (and possibly anomalous) Yang- 
Mills theory. It is rather obvious that the truly new element in these types of actions shows up when 
a physical axion-like particle is induced in the spectrum. In the absence of this, the bFF has dubious 
meaning, since this term does not cancel the anomaly, as emphasized by Preskill long ago [29]. Rather, 
it allows a better power-counting of the modified (anomalous) action. A justification of this point of 
view comes from the fact that an anomalous (and massive) Yang-Mills theory can be given a typical 
Stiickelberg form and a bFF interaction by a field-enlarging transformation 
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For this reason the only satisfactory potentials are those that either allow b to be part of the scalar 
sector (such as for V + V') or, alternatively, when they allow, under certain conditions that we are 
going to discuss next, a mixing between the CP-odd Higgs components and the Stiickelberg. 

With these motivations in mind, we move to the case of the new superpotential. 

7 Scalar mass terms, the scalar potential and the mass of the gauge 
bosons 

Let's now move to a discussion of the other sectors of the theory, starting from the scalar one. The 
Lagrangean for the scalar mass terms is given by 

CsMT = -MlDL-m%R^R-M^Q''Q-ml^u}iUR-ml^D\iDR-mlH^^^^ 
-mlHlH2 - mlS^S - {axSHi ■ H2 + h.c.) - {aeHi ■ LR + h.c.) 
-{aaHi ■ QDr + h.c.) - {a^H^ ■ QUr + h.c), (24) 

where Ml, MQ,rnR,rnijj^,m]yj^,mi,m2,rns are the mass parameters for the explicit supersymmetry 
breaking, while ae,a\,au, arc coefficients with mass dimension one. 

The computation of the Lagrangean containing the soft-breaking terms Lagrangean is, as usual, 
split into the scalar and gaugino mass terms 

C^Soft = C,SMT + C,GMT + (^bV'b + ^bV'b) , (25) 

where Mb is a mass parameter for the axino V'b- The gaugino mass terms given by 

C'GMT = —-^Mg (AgaApa -|- \ga\ga) — -^Myj (A^^iA^^i -|- A^iA^i) 

—My (Ay Ay + Ay Ay) - ^Mb {Xb^b + XbXb) , (26) 

where Xgo. , Xga are respectively the left- and right-handed gauginos of the SU (3) sector, X^yi , A^yi are 
the left- and right-handed gauginos of the SU{2) sector and Ay, Ay are the chiral gauginos of J7(l)y. 
The Mg,Mu),My,Mb mass terms are the SUSY breaking parameters for SU{3), SU{2)w, U{1)y 
and U{1)b respectively. Once we have imposed the equations of motion for the F-terms the on-shell 
Lagrangean is given by 

jCaux-F = -vMHiRR^ - ylHlH2URUl - vIhIh^DrD^j^ - \XHi ■ H2\^ 
-\XS\\hIH2 + hIHi) - ylDl^DRQ^fQ - yl^LRR} 
-ylURUlQ^'Q - Ay„ [sQ'^HrU]^ + h.c^j - Xya [sq''H2D\^ + /i.e.) 

Aye {sVH2R^ + h.c}j - yayu {uJ^hIh^Dr + h.c.) - y^ya {d^rQ^LR + h.c)j , 

(27) 
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where the coefficients ye,yu,yd come from the Yukawa couphngs of the superpotential, while the D 
terms are 

Cau.-D = -fiL^r'L + Q^t'Q + hIt'Hi + Hlr'H^f - ^{Q^T'^Q + Ur^T'^Ur + DR^T^DRf 



„2 

-^[BlL^L + BrR^R + BqQ^Q + BuU^Ur + BdD^^Dr 

+BH,HlHi + BH,HlH2 + BsS^Sf, (28) 

where Bl, Br are the charges of the leptons chiral superfields, -Bg, Bjj, Bd are the charges of the left 
and right chiral superfields of the quark sector and Bhi , , Bs are the charges of the two Higgs 
doublet and of the extra singlet respectively. 

7.1 The scalar potential 

The study of EWSB in the case of these models proceeds similarly to the USSM [20]. 
The scalar potential is given by 



V = \XH,-H2f + \XS\\m^ + \H2\^) + l{gl+g^){HlHi-HlH2f 

o 

2 2 

+ i^{BH,HlHi + BhMH2 + BsS^Sf + ^-\h\H2\'' + ml\Hi\'' + m^jFal 



2 

„2 I ol2 



+ml\SY + {axSHi- H2 + h.c.). (29) 
We introduce the following basis 

Hi = —\ ^ ^ \ , H2 = ^\ \ i \ , S = ^ Re5 + 1 Im5 , 

V2 \ ReH- + i IxnR- ) ^/2\ ReH^ + i ImH^ J V2 

(30) 

where in correspondence of the minimum value of the potential we use the following parametrization 
for the Higgs fields 

As usual, we require the existence of a stable vacuum imposing the conditions 

mfvi + ^X'^vi{v2 + vl) + -^axV2Vs - ^vi{vl - vj)g^ 

+LIBhMBhA + BhA + Bsvl) = 0, (32) 
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+\91Bh^V2{BhA + BhA + Bsvl) = 0, 



(33) 



-^axViV2 + msvs + ^X^vsv'^ + ^gBBsvsiBH2V2 + + Bsvs) = 0, (34) 

where again ax is a mass parameter of the model. 
7.2 Mass of the gauge bosons 

The Lagrangean that describes the contributions to the mass of the gauge bosons is given by 

1,.. . ,2 



Cq = \V^Hi\' + \V^H2\' + \V^S\' + - {d^lm b + MstB^) 



(35) 



and involves, beside the two higgses, the SM bosonic singlet of S, the bosonic component of the 
Stiickelberg axion, b, and the Stiickelberg mass Mgt- Collecting the quadratic terms we obtain the 
contributions to the gauge boson masses which are given by 



92 / ,2 



^3 92gY ,2 



+^(.? + vDA^^Al + '-^{BhA - BnAWl^' - '-^{BnA ' BhA^B^ 



+ ^-f{BlA + BlA + BhDB^B^^ + -M^tB^B^. 



(36) 



Using the interaction basis of the gauge field components (W^, A^, 5^) we obtain the corresponding 
mass matrix, which is given by 



'■gauge 



.92521^2 gy ^2 



9Y 



XB 



92 



XB 



9y 



XB 



+ 



Ml 



(37) 



where 



XB = 9BiviBH,-viBH,), NBB=9l{BjiA + BHA + Bhs)^ v'^ = vl + vl 



(38) 
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Performing the diagonalization we obtain the rotation matrix 



O 



susy 



gy 
g 



S2 

g 





gxB V2 



gy(/i~V/f+4x|g^) 



g2(fl-y/f!+4xlg^) 
\ gj2[Ag^xl+fl-h V/f+V^ gj2[4g2:^,2^+/2_/, ^7^^ 



qXB\/2 



ig'xl+n-h^iHig^xl] J 



(39) 



which acts on the interaction basis as 









Z 


^ susy 








\ B 1 



(40) 



and where we have defined g = gy + 92 fi — ^^st "5^^^ + ^bb- 

We obtain one nuh eigenvalue corresponding to the photon, while the masses of the physical Z 
and Z' are given by 

M| = i (^4M2 +g^v^ + Nbb - -g'^v^ + Nbb? + 



AMI +g^v^ + Nbb + J (4M|t -g^v^ + Nbb? + Ag^x% 



(41) 



Compared to the non-supersymmetric case |8j, the corrections to the masses of the gauge bosons 
involve also vs, which is implicitly contained in Nbb- 

7.3 The charged and the CP-even sectors of the scalar potential 

The description of the charged sector of the model is performed using the standard basis (Rei?;^ , Re-ffj" ) . 
We obtain the following mass matrix 



(42) 



The same mass matrix is obtained in the basis (— Im/f^, Im/7^ ). We have one zero eigenvalue corre- 
sponding to a charged Goldstone boson and a mass eigenvalue corresponding to the charged Higgs 



V2 Vl 



1 



1 



VS 



\ 75 ^1^2 - TT-^ ViV2 + . 



V2 



(43) 



17 



In the analysis of the CP-even sector we use the basis (ReiJ°, ReiJg, ReS). We obtain the matrix 
elements 



(A^L)l2 



(■A^ej33 



1 



(sl^li +9y + qI) vl - ax 



V2VS 

V2vi 



—Bh^Bh2 + a 1 viV2 + «a;^ 

ax^+I^^BH,Bs + \^^ v,vs 

1 / 2 d2 , 2 , 2\ 2 '"2VS 
1 {9bBhi +9y + 92) V2 - '^>^~^^^ 

+ [^Bh.Bs + A^^ V2VS 



-ax 



V\V2 



1 



with the other terms obtained by symmetry {Aii2 = A^2i)Ctc.). The matrix has in general three 
massive eigenvalues corresponding to the three neutral Higgs particles {H^,H2,H^)- 

7.4 The Neutral CP-odd sector and the cixion 

The key sector that is responsible for the presence of a physical axion is the CP-odd one. Choosing 
the basis given by the components (ImS", ImiifJ', ImiJg), our superpotential with an extra singlet gives 
the mixing matrix 



a\ 
V2 



V1V2 
VS 


V2 


Vl 


\ 


V2 


V2VS 
VI 


VS 




Vl 


VS 


VlVS 
V2 


J 



(44) 



Diagonalizing this mass matrix we can identify the orthogonal transformation 0°'^'^ from the interaction 
to the mass eigenstates which is given by 



/ Im5 \ 
Imi^o 



O 



,odd 



GO 



(45) 



A simple analysis gives two null eigenvalues, corresponding to two neutral goldstone bosons, and one 
physical state, which is identified with a massive neutral Higgs boson 



m 



2 _ ax f nV2 ^ VlVS ^ VSV2 
V2\VS V2 Vl 



Hi 



(46) 
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Prom the diagonalization procedure we obtain 



O 



and the states are given by 



odd 







vs 



VlV2 







V2VS 



V-LVS 



(47) 



vilmH^ — vslmS 



^4^ = 



(48) 



where and G2 arc two Goldstone modes, while is the physical Higgs. 

Having identified the goldstones of the potential in the CP-odd sector, the parallel identification 
of the goldstones of the massive gauge bosons after EWSB is performed by an analysis of the bilinear 
mixings. In fact, from the Lagrangean density we can extract the following derivative coupling terms 



£^DC 



-g2W^d''GY - -gvAl d^^Gy + -gsB^d^^GB 
where we have defined 

Gy = {vilm Hi - valm H^) 

Gb = (BHiVilm + Bh2V2^ + fig-y^Im S) + Im h 



(49) 



9B 



(50) 



which can be rotated onto the basis (^/l, Z^, Z'^) using the Of^^gy matrix 



B, = O^^Al + OizZ, + 0^z,Z'^ 
to obtain the expression for £dc in terms of physical states 

Cdc = MzZf.d^Gz + Mz'Z'^d^Gz'. 
The two goldstone modes corresponding to the physical massive gauge bosons are given by 



(51) 



(52) 



MzGz = -A 



2XB 



h + sJft + ^g^x%\-v^gBBH, 
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+A 



V2 
2XB 



fi + \/fi + ) + V2gBBH., 



ImHl 







+BsgBVs Aim S + 2MstAlmb 



(53) 



Mz'Gz' = A' 



2XB 



ft + ^9^xl-h]+vigBBH, 



-A' 



V2 

2^ 



/^TV^ - /l ) - V2gBBH2 



+BsgBVs A!lm.S + 2Mst A' Im b 
where we have defined the following coeflacients 



A 



fi 



A' 



1 

t: + 



(54) 



(55) 



It is simple to observe that Gz and Gz' arc orthonormal. At this point, a simple counting of the 
physical degrees of freedom before and after EWSB can give us a hint on the properties of this model. 

Before EWSB we have ten degrees of freedom: two for A^j^, two for W^, three for S^, two for the 
Higgs fields ImH^ and ImiJ^ and one for the singlet Im S. After the breaking, we are left with two 
polarization states for the physical photon, three degrees of freedom for the Z and the Z' respectively, 
one neutral Higgs state H2 and one physical state which we are going to identify as the axi-Higgs. 
Therefore we can build this new physical state requiring its orthogonality with respect to the basis 
{i?4, Gz, Gz'} where H^, identified as the physical direction of the potential, clearly belongs to the 
CP-odd sector. We start from the following linear combination 



X = 6ilm + 62lm + b^lm S + b^lm b 
and we determine the coefficients 6i, . . . , 64 by the following constraints 

li = 63^1^2 + b2Vivs + biV2Vs = 0, 

Y2 = Ab^MstXB + 2b-iBsvsgBXB - bivi{fi - 2BH,gBXB + yj fi + Vx|) 

+b2V2{fl + 2BH,gBXB + ^/f + VxD = 
Y3 = Ab^MstXB + 2bzBsVsgBXB + b2V2{fl + 2BH2gBXB - 



which give 



(56) 



+6l^^l(-/l + '2BH,gBXB + Jff + Vx|) = 0, 



(57) 



bi = 64 



2M, 



St 



V1V2 



gBBs {vlvl + v'^v^g) 
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b2 
b3 



M. 



St 



VIV2 



M. 



St 



where the coefRcient 64 is constrained by the normalization of the eigenstates. The physical axion will 
be given by 



X 



[2MstVivl ImH^ + 2Mstvjv2 ImH^ - 2Mstv'^vs ImS + Bs gsiv^vl + vjvl)lm b] 



(58) 



where the new identified state has a nonvanishing projection over the Stiickelberg field. Re-expressing 
Im6 in terms of x and the goldstone modes of the massive gauge bosons, we discover that the axion- 
like interactions (Wess-Zumino terms) mediated by the Stiickelberg field can be rotated over Xj giving 
trilinear vertices of the form x^i ^ where I and J denote the physical gauge bosons. 

The rotation matrix O^usy that rotates the physical components and the goldstones in the CP-odd 
sector takes the form 



Gz 
G' 



susy J 



\ X J 

where all the entries are defined in Appendix B. 



/ Im \ 
Im 
Im 5* 
Im b 



(59) 



7.4.1 The Bs = case: no physical axions 

In the case Bs = 0, corresponding to a singlet of the entire gauge symmetry, we can proceed in 
the same way, obtaining, however, a different result compared to the previous case. In this case the 
general structure of the scalar potential can be modified by introducing linear or cubic terms in S, 
corresponding to the same structure of the nMSSM or of the NMSSM, with an additional U{1)b 
symmetry. Adding a linear term we obtain [§ 



V 



IXHi ■ H2 + 



m 



12 |2 



A 



+ \XS\'{\Hi\' + \H2\') + -{gi+g{r){HlH_ 



B\B2f 



+^BI^{hIh, - HlH2f + %h\H2\'' + m\\Hi\'' + ml\H2 

o Z 

+m||S|2 + [axSHi ■ H2 + tgS + h.c), 



(60) 



■^At this stage we do not consider a cubic term in S in order to avoid the problem related to the formation of 
cosmological domain walls (see [22], [31], [32]), though even in this case one has two Higgs bosons and one Goldstone 
mode in the CP-odd sector. 
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where we have introduced the mass parameter m^g/-^ " which is the coefficient of S in the nMSSM 
superpotential - and ts , which is the coefficient of 5 in the soft breaking Lagrangean and has mass 
dimension three. Notice that we have used the condition Bh-^ = —Bh^. In this case, in the basis 
{ImS", Imff^jlmi^g}) the CP-odd mass matrix is given by 



odd 



(61) 



This sector provides two physical Higgs states and one goldstone mode of the form 



i 



GImssm = J-^ (film - "2lm . (62) 

The other goldstone mode is obtained from the derivative coupling of the Stiickelberg term {B^dfjlm b). 

Thus, from the derivative couplings, once we have performed a rotation on the physical basis, we 
obtain the two orthogonal Goldstone modes Gz, Gz' corresponding to the Z and the Z' bosons, which 
are a linear combination of Im b and of the Goldstone mode obtained from the CP-odd sector, 

Gz = aiG^MSSM + a2lm b, Gz' = aiG^MSSM + "2!™ ^ (63) 

where the coefficients ai . . . ,a'2 are not given in an explicit form for simplicity. 

In this case the number of degrees of freedom before the symmetry breaking is again equal to ten. 
In fact we have two for three for B, two for Y and finally Im H^, Im H2 and Im b. After EWSB we 
are left with three degrees of freedom for the Z, three for the Z' , two for the photon and two neutral 
higgs states, which are physical. Therefore we do not have Higgs-axion mixing. 



8 The sfermion sector 

Coming to the scalar fermion sector (sfermions), the Lagrangean in terms of component fields is given 

by 

CffSjM ^ -\y^[S^HlLR + SVH2R^] - A yd[S^HlQDR + SQ^H2D\^} 

-XVuIS^hIqUr + SQ^HiUj^] - yl[HlHi{VL + W R) - h\l{hIl)^] 
-yj[HlHi{Q^Q + D^^Dr) - hIq{h\Q)^] 
-yl[HlH2{Q^Q + uIUr) - hIq{hIq)^] - MIDl - mlR^R 
-MlQ^Q - mfj^ulUR - ml^D^DR - {aeHi ■ LR + h.c.) 
-{adHi ■ QDr + h.c.) - {auH2 ■ QUr + h.c.) 
^The same goldstone mode can be obtained from the NMSSM scalar potential [33) . 
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-^{Dl - 2&R - igtg + ^ulijj^ _ Id^Dr + hIh, - HlH^f. 



(64) 



In the presence of an extra U{1)q an additional piece coming from the D-terms must be added to the 
sfermion Lagrangean and it is given by 

2 

'^'tfer'' = -y^BlL^L + Br&R + BqQ^Q + BuUI^Ur + BdD^rDr 

+BhMHi + BhMH2 + BsS^Sf. (65) 
After spontaneous symmetry breaking we get 

4/er = -l\vsyeV2[L^R + L^^ R^] - ^XvsydV2[Q^ Dr + Q'^-DJj] - ^XvsyuVilQ^UR + Q^^UU 

-^yy,[P^P + R^R] - ^yjvjlQ^^Q' + D^^Dr] - \ylvl[Q'^Q^ + uIUr] 
-MlVi - mlR^R - MlQ^Q - ml^ulUR - ml^D^^DR 
-{ae^L^R + h.c.) - {aa^Q^DR + h.c.) + {au^QiUR + h.c.) 

o 

- v'2){L''L - 2R^R - igtg + ^ ~ ^ _ ^t^^^^ 

„2 

{Bh.vI + BhA + Bsvl) {BlL^L + BrR)R + BqQ^Q + BuUIUr + BdDIDr); 

(66) 

here and in what follows superscripts on L and Q specify the doublet components. 
In the basis (Z^, w), the entries of the mass matrix are given by 



(^L^i?)ii = yll^l + ligl - ghivl - + ^Bl{BhA + BhA + Bsvl), 

1 Vl 

{Mi2^r)i2 = {Mi2^r)21 = -XvsyeV2 + ae-j=, 



(M^, b)22 = lyX + ml- -gl{vl - vl) + ^BRiBn^ + BhA + Bsvl). (67) 



The former matrix can be diagonalized through a rotation defined by 

tan2^^,^ = {XvsyeV2 + aeV2v,) _ ^^^^ 

' rnl-Ml + \{gl - 2,gl){vl - vl) + '-f{BR- Bl){Bh,vI + Bh^vI + Bsvl) 
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The eigenvalues have very lengthy expressions and we will omit them for brevity. The three eigenstates 
are given by 

h = cos 6^2 ffL^ + sin 6 12 fjR' 
h = — sin 9i2ffL?' + cos Oi2fiR^ 

h = L\ (69) 

The mass of is given by 

Ml = + ghivl - vl) + ^^Bl{BhA + Bh2vI + Bsvl). (70) 
Using the two basis (Q^, -D]j) and {Q^,tj\^, the mass sector of the squarks can be written as 

Cs,uark = -(g^t ) M^.^^^ ( 1^ J - ( Qlt l}R)M^,^^i^\)^, (71) 



where the Mq2 matrix is defined as 



(Mq^^^Jh = ]^ylvl + Ml- ^-{gl + \gl){vl - vl) + ^-^Bq^BhA + Bh2vI + Bsvl), 

1 V\ 

(^Q2,i5«)22 = \yyi + ml^ - Y2^Y{yl - vl) + ^Bd^{Bh,vI + Bh2vI + Bsvl), 
while for the Mgi matrix we get 

= ^ylvl + + lial - prM - vj) + ^Bq{BhA + Bh2vI + Bsvl), 

1 V2 

y^Q^ijRn^ = (A^Q^^/«)2l = i^>^vsyuvi - 

(Mqi j>^)22 = \ylvl + ml^ + ^gl{vl - vl) + ^Bur{BhA + Bh2vI + Bsvl). (72) 
The Mq2 matrix can be diagonalized using 

qi = cos0Q2_^^Q2 + sin0Q2,^^£'t 
q2 = -sin0Q2^^^Q2 + cos0Q2_^^£>lt, 

where the 6q2 jj^ angle is defined by 

tan2^.,^ = {XvsyaV2 + a,V2v^) _ ^^3^ 

' " ml^-Mi^ + ligl - yi)ivl - vl) + f (B^,, - BQ)(5^,t;f + Bh2vI + Bsvl) 
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Again, we omit the explicit expression of the eigenvalues since they are quite lengthy. The Mqi 
matrix can be diagonalized by the following choice 

q4 = ~ sin ^^^Q^ + cose Q.^^^UI, 

where 6qi is defined by 

tan 20.,^ = {XvsyuV2v^- a^V2v,) ^^^^ 

^ ' " - - ligl - lgl){vl - vl) + '-f{Bu, - Bq){BhA + BhA + Bsv^) 

Using the parameter values specified in the numerical analysis of the neutralino sector, typical values 
for sfermion masses are around a few TeV. 



9 Wess-Zumino count erterms and Chern- Simons interactions 

The cancellation of the gauge anomalies in these supcrsymmetric models are obtained by the introduc- 
tion of axion counterterms. The supersymmetric form of the corresponding Lagrangean introduces, 
beside the usual bosonic contributions of the form bF A F additional interactions between the axion 
and the gauginos and between the axino, the gauge fields and the corresponding gauginos. It is given 

by 

Cc = - J d^e!^ ^bGTr{gg)h + ^bwTriWW)h 

+bYhW^W^^'^ + ftfibW^W^^'" + bYBhW^W^''' 6(6^) + /i.e.} . 



(75) 



whose general e Expanding this expression in component fields using the WZ gauge we obtain 
jOc = -lbG ei^'^P'^Gl^Gl, Im6 -^-bw e^'^^- W^, W;, Im6 
-lbye^-P'^F]:,Fj,lu.b - \bBe'^^P'^ F^^F^^l^b - ^byBe''''''' F^^F^^lmb 

+ bG[Imb^{Xgaa''D^Xga) - _!^Vb^(Apa(7'^a'^G« J + ^F^^{XgaXga) 

+-Lv;b^(A,.DS) + h.c] + bwilmb^iX'^a'^D.X'^) - -l^^^liX^^a^a^W^^^) 

+\F^\i^wi^wi) + ^V^b^lAvyiD^) + /i.e.] + bY\\T^lbXYa^D^,XY - ^VbAya'^a'^Fj^ 

11 - i 

+-FbAyAy + ;^V'bAy -Dy + /i.e.] + ^^[Imb Asa'^D^As - ^VbAsa'^a^i^^ 

+^FbABAB + ^VbAs Db + /i.e.] + 6yB[(Im6Ay(7^a^AB + ^FbAyA^ 
z y 2 z 
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+ -^V'bAy Db - -^Xyo^a^'F^M + (y ^ 5) + h.c\. 



(76) 



where we have additional contributions for the cancellation of the C/(l)B5C/(3)S'f/(3) anomaly, which 
are typical of this model and are not present in previous similar formulations [lOj . 



9.1 The Chern-Simons Lagrangean 

As we have mentioned above, the Chern-Simons Lagrangean describes the freedom to re-distribute 
the anomaly in the trilinear gauge interactions of AW and AAA type. In a bottom-up description 
of these models this freedom is equivalently formulated in terms of external Ward identities on the 
anomalous vertices. The corresponding Lagrangean is similar to the one given in [TO], now with the 
addition of the gluonic terms. It takes the form 

Ccs = - j d"^^ {ci {YD'^B - BD'^Y)W^ + h.c. 
{YD'^B - BD''Y)Wl + h.c. 



-C2 

-C3TV 
-C4Tr 



{WD'^B - BD''W)Wa + -WD'^BD^lDaW, W] + h.c. 

(GD^'B - BD^'G^ga + -GD'^BD^lDaG, G] + h.c. 

6 



(77) 



where the coefficients ci . . . C4 will be determined by the generalized Ward identities of the model. 
Expanding this expression in terms of component fields we get 



C 



cs 



-c^e^'^P'^B^Y.F^^ + C2e^"'"^S^y.F/„ + cge'^^^'^B^TV { W^F^^ - -W,[Wp, W„] 



+c^e''''P''B^Tr [G,Gp^ - -G,[Gp, G^]j - ci (Ab^^AbA^ - Xbct^^XyB^ + h.c.) 

+C2 (Ayfj^AyS^ - AyfT^Afi^^ + h.c.) + C3 TviXwa^" XwBf, - Xwct^'XbWp + h.c.) 
+ C4 TiiXga%B^, - Xga^XeGf, + h.c). 



(78) 



The role of the Lagrangean is to redistribute the anomaly among the three anomalous vertices when 
the symmetry of the interaction is not enough to fix the partial contributions to the anomaly uniquely. 



10 Generalized broken Ward identities 

The anomaly cancellation mechanism for this supersymmetric model proceeds as in [51 [30} [Ml [36} [371 
I38j . where a detailed description of some physical cases can be found. The resulting anomalies must 
be cancelled in the abelian sector BBB, BYY, YBB and in the non-abelian SU{2) and SU{3) sectors. 
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If we start by using a parametrization of the one-loop trilinear gauge interactions with a symmetric 
distribution of the AAA anomaly vertex (Aaaa), in which we denote with —ks = ki + k2 the incoming 
momentum with the A index and with ki , k2 the outgoing momenta, with indices /x and v respectively, 
we can introduce generalized Ward identities in the momentum space as defining conditions on the 
model. We obtain 



k3,xABBBA\^XAik3, kiM) - \hB e'"""^ki,ak2,p - 2m/ A^'^^ = 0, 



(79) 



for the BBB case, and analogous conditions in the other sectors. The expressions of Aaaa, ^bb and 
similar are given below; ruf denotes the mass of the fermion in the anomaly loop. 

Other two Ward identities are obtained by a cyclic permutation of the momenta. Also, notice that 
in this specific case we do not have Chern-Simons interactions in the defining condition. For a BYY 
triangle we have 



k 



'3,A 



ABYYA\^XAik3, ki, k2) - C2e^^''"{ki 



k 



■2)a 



1 



bYe'"""^ki,ak2,p - 2m/A^V = 0, 



^1,^ [ABYYAXXAik3,ki, k2) - C2£^''""(fel - k2)a\ - 2mfA^^Y = 0, 

k2,. [ABYY^AAAiks, k^M) ' C2e^>""'{ki - k2)c] - 2m/A^'^^ = 0, 
where the tensor structure of the triangles is given below. For a YBB triangle we have 



(80) 



k3,X 
k2,u 



^yBBAX(fc3, ki, k2) - cie^^'^"(fci - fca)^] - 2m/A^'^5 = 0, 

.4yBBA^Xi(fc3, ^2) - ci£^'''^"(fci - k2)a] - \hYB e'"""^k2,ak3,p - ^rrifA^B^B = 0, 

AYBBA^X^iks, kiM) - cie^^'^'iki - k 



2 a 



(81) 



where the coefficients ci, C2 are fixed by the BRST invariance under C/(l)y. The explicit form of the 
tensors A^'^^ and A^'^, in terms of Fcynman integrals, are given by 



^AAAi-^f ^^) = ^ [ dx 
^ Jo 



e[ki,X,fj,,i'] 



+£[A;2, A,/i, I/] 



l-x 



A{mf) 



dy 



A(m/) 



{ 



3 

A{mf) — m?j 



— ^ + k2 ■ k2y{y - 1) - xyki ■ k2 



ki ■ kix{x — 1) + xyki ■ k2 



+e[ki, k2, A, i']{k'^x{x — 1) — xyki^) 
+£[fei,fc2,A,/n](/c2y(l -y) + xyfci)} , 



and 



37r2 



2/3 



1 rl-'x 



^0 



dxdy 



A{mf) ' 



(82) 
(83) 
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where A(m/) = [mj + {y — l)?/^! + — l)xkf — 2xyki • k2\ ■ For Ayy and Ay^ we obtain similar 
expressions. The same relations can be reformulated in the mass eigenstate basis in terms of the 
physical gauge bosons Z and Z' . The structure of the (generalized) Ward identity in this case is 
shown in Fig. [H written in configuration space, where the first term corresponds to the anomaly, the 
second is the axion counterterm projected out on the goldstone Gz, and the third diagram describes 
the mass corrections due to the coupling of the goldstone to the massive fermion in the loop. In the 
chiral limit, obviously, the third term is absent. 




a) b) c) 



Figure 1: The generalized Ward identity for the Zjj vertex in our anomalous model away from the chiral 
limit. The analogous STI for the SM case consists of only diagrams a) and c). 

The generalized Ward identities for the case U{1)b SU (2) SU{2) have similar expressions, while the 
case U{1) B SU{3) SU (3) requires a further comment. As a matter of fact, in this case the higgsinos do 
not circulate in the loop, but the BGG triangle exhibits an anomaly when Bs 7^ 0, (see Eq.Q). For 
the same reason we do not have a BGG anomaly in the MLSOM [5] (Minimal Low Scale Orientifold 
Model) case when the Higgs charges under U{1)b are equal. 



11 Z decay into four fermions: Chern-Simons interactions 

One interesting signature of trilinear anomalous vertices involving three anomalous gauge bosons can 
be investigated in the decay process of the Z/Z' into four fermions by the mediation of two extra 
anomalous currents. This kind of process is phenomenologically relevant since it is sensitive to the 
presence of (at least) two or more extra anomalous U{1). As a matter of fact, in the MLSOM 
(non supersymmetric case) in the presence of an abelian symmetry given by Gi = U{1)y x U{\)b 
where B is anomalous, the off-shell effective vertex does not contain any Chern-Simons interaction 
by construction. If we take, for instance, the triangle {ZZ' Z')^ some of the relevant effective vertices 
coming from the interaction eigenstate basis which have an anomalous component are (BBB) and 
(YBB). In the BBB case the Chern-Simons interaction vanishes trivially, while in the YBB case 
the corresponding Chern-Simons counterterm must be "absorbed" in a redefinition of the triangle in 
order to ensure the BRST invariance. Equivalently, the YBB vertex does not allow a partial anomaly 
on the Y leg, since there is no axion for Y . An analysis of the anomalous trilinear interactions in the 
context of the MLSOM can be found in |36| . 

In the presence of multiple anomalous C/(l)'s (such as U{1)y x U{1)b x U{1)bi) the situation 
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is quite different. The Z decay into four fermions can be mediated by two different extra neutral 
currents and the off-sheh vertex can be of the type (ZZ'Z"), while from the interaction eigenstate 
basis a contribution BB'B' appears. A simple inspection of the gauge invariance of this vertex shows 
that a Chern-Simons interaction can not be absorbed into a redefinition of the BB'B' triangle. 

A symmetric distribution of the anomaly on the BB'B' triangle, with outg oing momenta h-^ , h2 
and incoming momentum k, fixes the Rosenberg parametrization as follows @ 

"^AAA = (-A^i • ^2 - Aekl - -^^/^j'^] + i-^Aki ■ k2 - A^kf + y)^[^2, A, 

+A3k'^e[ki,k2,X,u] + A4kt^e[ki,k2,X,i'] + A5k'[e[ki,k2, X, fJ'] + AQk^£[ki,k2, X, fj], (84) 
thus, we have a partial anomaly equal to ^ on each Lorentz index 

K'^AAA = ye[^i> ^2, A, y] 

klT\^'AA = -^^[kiM,\lA- (85) 

The generalized Chern-Simons interaction allowed by the presence of multiple anomalous C^(l)s can 
be formally written as 

V^^s" = ^, - ^2 ) + 4'^^[A, M, a]{k^ - k^) + a(f)e[A, ^, a]{k^ - kf) (86) 

(i) 

where k^ = —k and the coefficients i = l7 2,3 depend on the model and satisfy the relation 
On ^ + ttn^ + = a„. Therefore, in the definition of the effective vertex the contributions coming 
from the Chern-Simons interactions appear explicitly and spoil the symmetric distribution of the 
anomaly on BB'B' . Moreover, the cancellation of the anomaly is ensured by the presence of the 
WZ interactions, which are constrained by the BRST invariance of the model. For example, the 
computation of the diagrams described in Figs. [2] and [3] gives 



f = e\k) m^AA + yc^" 



a""" - -TTT- 1 ,2 :.2 ^(g3)r.^^(g4) 



(87) 



Mj„ k^- M|„ 

where we have indicated with T^i the generic Lorentz structure of the fermion coupling to the extra 
Z' /Z" . For instance, the Chern-Simons contribution gives 



Tcs = a(^)e[A, ^, u, ki - ^2] + a^'^h[X, fi, u, k2 - k^] + a^-^h[X, fj., k^ - ki] 



X 



-1 



"'^■'r''"'''^'"fe'r''"fa' (t;-M|,)(>.^-M|„) - 

The detection of these interactions is rather difficult experimentally, given the low production rates 
due to the large mass of the extra Z', currently bound to be larger than 900 GeV. 



We have defined a„ — and we use the notation e[fci, ^2, ^, v\ = e'^^^^k\^oM\ 
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(a) 



(b) 



(c) 



Figure 2: Redefinition of tlie effective trilinear vertex including the Chern-Simons interactions. 






Figure 3: Decay of the Z boson into 4 fermions plus the Chern-Simons contribution. 



12 The Neutralino sector 



Moving to the neutralino sector, here the mass matrix is 7-dimensional because of the presence of the 
axino, the singlino and the B-ino in the spectrum. In the Bs ^ case we obtain 



1 11 IVi ~ IV2 ~ 2 I'Vl ~ I 

--MyjX^zXwi — -My Ay Ay - -MbXb^B + -^92^W^H-^ - -j=g2\w3H2 — -^gYXyHi 

+ '^gY\YHl + '^gBBH.XBHl + '^gBBH^XsHi + '-^gBBsXB~S - X vsHIhI 



-A viSHl - X V2SHI + ^i^h^B - 2^bV'b^b + h.c, 



(89) 



where My^ , My , Mb , Mb are mass parameters and the term Xvs/ \/2 plays the role of the ;u-term; notice 
that A is a dimensionless parameter. We have indicated with A^3, Ay, Xb the gauginos of W^, , B 
respectively and with the SUSY particle associated to b. The fields HI and (i = 1,2) denote 
the supersymmetric partners of the two Higgs doublets, while S is the SUSY partner of the extra 
singlet S. 
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In the basis (— iA^a, —iXy, —i^B, -f^i; -f^l ' "^V'b) the mass matrix takes the form 











-f52 


f52 











My 







-fffy 














Mb 




-^gBBH2 


-^gsBs 


V2 




f5y 









"^73 



























-"fgBBs 




X ''I 

"^73 














Mst 











Mb 



(90) 



V 

that will be analyzed numerically in a section below. 
12.1 A preliminary choice 

A preliminary choice [20j which allows to simplify the structure of the 7x7 neutralino matrix is 
made by setting M^ = My = Mb = Mb = A = 0. In these conditions the diagonalization is rather 
straightforward and we obtain three null eigenvalues. The first corresponds to a physical pure-photino 
which is obtained from the rotation 



A^ = sin Ow^w^ + cos^vK'^yj 
^ZsM ~ cosOw^w^ ~ singly Ay, 



(91) 



where Xzsm is intermediate unphysical state. The second state, corresponding to a null eigenvalue, 
is given by a mixture of Higgsino and axino states 



X'2 



M. 



St 



M, 



St 



-Hi + V'b, 



(92) 



2gBViBs 2gBV2Bs 

while the third is a pure Higgsino state which corresponds to the SUSY partner of and it is given 
by the expression 



Vi V2 



(93) 



The other states corresponding to the non-zero eigenvalues are complicated combinations of higgsinos, 
gauginos (A^g^^, A^) and the axino. 

Notice that in our treatment we are considering for simplicity a real-valued neutralino matrix. 
In the most general cases - for example in some CP-noninvariant theories - these matrix elements 
are complex and they may contain phase factors which are physical and can not be eliminated by a 
redefinition of the fields. 
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z z z' 

(a) (6) (c) 



Figure 4: Trilinear interactions between x ^-nd the neutral currents 



13 Supersymmetric interactions of the axion with the neutrahnos 

In this section we proceed with a study of the basic tree-level interaction vertices involving the phys- 
ical axion (axi-Higgs). Analyzing each sector of the whole Lagrangean we have different types of 
interactions involving the axi-Higgs. 

First of all, from the counterterm Lagrangean we have trilinear interactions obtained by rotating 
the WZ counterterms on the physical basis, which formally give terms of the type 

£75 ^livp^ ryahel ryahel , I r> ^iivpfy ryiahel ry label , , n ^livpfy ryahel ly label , ((\A\ 
xZZ = tiie>^ ^ Z^^ x + Z^^ X + -Ks e'^ Z^^ X, (,y4j 

where for simplicity we have indicated with Ri,R2,R3 the coefficients which appear in front of each 
vertex. These include the rotation matrices, the coupling constants of the gauge groups and the 
coefficients coming from the anomaly cancellation procedure. We omit their explicit expressions since 
they are not relevant for this discussion. Notice that in this case only the abelian part of field strengths 
contribute to the counterterms for the neutral currents and that Z^'^'^^ = df^Zy — dyZ^. The interactions 
coming from these terms are shown in FigjH 

From the axion Lag rangean C^axion obtain quadrilinear interactions between X; the neutrali- 
nos/gluinos/charginos, the neutral/charged gauge bosons and trilinear derivative interactions, illus- 
trated in FiglHMl In fact, by a careful inspection of Caxion we find 

^^fon""'' = R^x 'xrx"^ z, + R^x 'arc G. + R'^x ~x^r^x^ wj + {z^ z'] , (95) 

while the derivative trilinear interactions are given by 

^ifon = X xlT^d^x] + R''^'^ X &l^d^G + R^^ X X^T^d^x"^ , (96) 

where P'^ indicates that we can have vector or axial- vector interactions. Trilinear interactions between 
one neutral current and two axion-like particles can be obtained from £quad and have the form 

^QvIb = R""""^^ X 9^ H^Zp + /?xH±14^T ^ fj±^T + ^ z'}; (97) 
to these terms correspond the interactions shown in FiglT} Analogously, the quadrilinear interactions 
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Figure 5: Quadrilinear interactions involving x, charginos/gluinos/neutralinos and a gauge boson. 




(a) (6) (c) 



Figure 6: Derivative trilinear interactions between x and charginos/gluinos/neutralinos. 

between two axion like particles and two neutral gauge bosons are given by (see Fig. [8]) 
-Cquad = Ri^ xxz^z^ + Rl^ X H^Z^Z'^ + i?f ^' X X z^z'^ + flf ^' X ulz^z'^ + {Z^ z'} 

(98) 

where, again, we have introduced the coefficients Rf^,R^^ containing the rotation matrices and the 
couplings, for simplicity. 

From the Lagrangean of the scalar mass terms Csmt we obtain the following trilinear interactions 
involving the axi-Higgs, the Higgs bosons coming from the scalar sector (CP-even, CP-odd, charged) 
and the sfermions 

CYmT""" = R"""' x' + R""' X HlH^ + R^^ xH^H^ + R'^^X If, (99) 

where Hf with i = 1, ... 3 indicates the physical Higgs states coming from the CP-even sector (see 
Figl9]). We denote with Cy^ the on-shell Lagrangean coming from the superpotential, once that the 
F-terms have been removed, and containing all the Yukawa-type interactions 

£yv; = Cyuk + + ^Yuk~F (100) 

where Cyuk represents the Yukawa interactions that do not contain the extra singlet S and are linear 
in ye,yu,yd, while £5 indicates all the Yukawa interactions containing S. Finally, with Cyuk-F we 
indicate those interactions that are quadratic in ye,yu,yd and in A. Then we have 

>Cy„fc = yee'^[-H\VR - HiVR) - H\VR - H'^VR - RH\V - RHi V'^] 
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X X 




(a) (b) 

Figure 7: Trilinear interactions between x, an Higgs boson and an electroweak gauge boson. 




Figure 8: Quadrilinear interaction involving x, two electroweak neutral gauge bosons and the CP-odd Higgs. 



+yue'^[-H\Q'tjR - AIq^uI - HIQWr - h'^Q'Ur - UrH\Q^ - UrAIq% (101) 
= -\ye[S^HlLR + SVH2R^] - \yd[S^HlQDR + SQ^/Za^J?] 

-\XS\\hIH2 + hIHi) (102) 



and finally 



Cvuk-F = -\XHi ■ H2\'^ - yl[L^WR + H\Hi{VL + R) 
-H\l{hIl)^] - yj[Q^QD\^DR + h\Hi{Q^Q + D^^Dr) - h\Q{h\Q)^] 
yj^UR + HlH2{Q^Q + ull 



-yl[Q^QUj,UR + HlH2iQ^Q + uj,UR) - hIq{hIq)^] . 



(103) 



From the Yukawa mass terms contained in Cyuk and in Cs we can isolate the pseudoscalar coupling 
of the axi-Higgs to the fermions and a quadrilinear scalar interaction with the sfermions 

^^uk-s = Rrii i^n'i^f X + Rt"x X If + Rf' "x Hi ff (104) 

where we have indicated with ipf the generic fermion and with / the generic sfermion (see Fig llOp . 
Quadrilinear axionic self interactions can be obtained from Cs and from Cyuk-F 

c^^'^ = R'<~\' + R'^'x'iil + R^'^x^H^H^ + R^Wulf + R^x^alf 

+R^^'^ X^H^H^ + R"""''^ XH^H^H^ + i?^^"^ xH^H^H^ (105) 
and are listed in Figllll 
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(a) (b) (c) 



Figure 9: Trilinear interaction involving x a-nd Higgs bosons/sfermions. 




(a) (6) 



Figure 10: Interactions obtained from >Cy^j,. 

14 Numerical Analysis 

In this section we present a numerical analysis of the neutralino sector. We have performed the nu- 
merical diagonalization of the 7x7 neutralino matrix and we have studied the eigenvalues dependence 
with respect to the free parameters of the model. Furthermore, since in this model the neutralino sec- 
tor exhibits an axino component due to the presence of Stiickelberg interactions, we have investigated, 
in the case of the lightest neutralino state, its mixing with the other states. In Tab. Owe have listed 
all the values of the parameters that we have used in our analysis. In our analysis we have followed, 
in spirit, the approach of Kalinowski and collaborators in [39]. In their paper the authors, who deal 
with the USSM, present two scenarios: in the first one they assume unified values for the gaugino 
mass terms and in a second scenario they consider with different values (arbitrary values). We refer 
to their analysis for further justifications and motivations of this choice. We have chosen tan/3 ~ 40 
and we have constrained the value of vi in order to be consistent with the value of the mass of the Zq 
boson, while the value of the coupling constant qb is 0.65. 

The values A < 0.7 and vs around 1 TeV are consistent with the MSSM value of the Higgs masses. 

The charges Bh^ and are free parameters because we have only four equations coming from 
the gauge invariance of the superpotential and eight charges to be constrained. One possible choice is 
Bh^ = -3/(2^/10) and Bh^ = -l/(\/TO), which is obtained from the EeSSM model [39]. 

In Figs ll2lfT5t we plot on the left-hand side the numerical value of the neutralino masses obtained 
from the diagonalization procedure as a function of the mass parameters Mgt, Mb, M]j, My , Mw and 
of gB and tan /?. On the right-hand side we plot the squared value of each component of the lightest 
neutralino state in order to establish which component is dominant, since every neutralino state 
appears as a mixture of the axino, the singlino etc. We can formally decompose the generic i-th 
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(a) (6) (c) 

Figure 11: Quadrilinear interactions involving x and CP-odd/CP-even/charged Higgs. 



neutralino state (i = 1, . . . , 7) in the basis {—iXw^, —i^Y, —i^B, Hi, S, —iiph} 

Xi = «ii ^Ws + «i2 Ay + Ojs Xb + an Hi + ai^ ij| + ate S + aa ip^ (106) 

and in the figures we indicate the square of each component I , where the lightest state 

corresponds to the i = 1 choice. From the left panel of Figsll2l and 1131 we observe that the value of 
the mass of the lightest neutralino state that is consistent with the current experimental bounds [35] 
is obtained approximately by varying the values of Mgt in the interval 1.7 2.5 TeV, while Mb and 
Mb in the interval 1-^2 TeV. In the right panel of Figs]12landll3lit is interesting to observe that for 
these values of the soft breaking parameters we have a tiny region beyond 1 TeV in which the axino 
and the B-ino components are almost coincident, the two higgsinos are dominant, while the singlino 
is the most suppressed component. For values of M^t, Mb,M]j below 1 TeV and beyond 2.5 TeV, the 
lightest neutralino is "mostly" singlino, while the M^-ino and the Y-ino components are suppressed 
and the eigenvalues appear to be non-degenerate apart from the states X2 ~ Xs- From the left-hand 
side of Fig. [T3] it is evident that all the eigenvalues do not exhibit substantial variations with respect 
to My,M^ and the heaviest states are non degenerate. In both cases (see Fig. [H] (b), (d)), the 
singlino component is the leading one. A similar feature can be found in the USSM case [39], where 
the singlino is always dominant with respect to the other components. 

Finally, in Fig. [15] we have analyzed the dependence upon the coupling constant gB, tan j3 and vs- 
In the left-hand side (a) the mass value of the lightest state starts to be greater than 50 GeV once 
qb > 0.4 and it is almost degenerate with X2- 

From the analysis of each component in the right panel (b), for gs less than 0.5 the main contribu- 
tion comes from the singlino, while the axino and the B-ino are almost degenerate and subdominant 
with respect to the H2 contribution. When gB becomes greater than 0.5 we have an inversion: the 
two Higgsinos are dominant and almost equal, while the singlino is subleading and the combination 
axino- -B-ino is more suppressed. 

As a consequence of our constraint on the vev vi , the eigenvalues dependence on tan (3 is weak 
(see Fig. [TJ] (c)), while we have a strong impact of low values of tan /? on the axino, B-ino and on the 
singlino components. Even in this case, with the choice of the parameters that we have made in Tab. 
[21 we can identify a small region in which the contribution of the singlino is highly suppressed. 
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In the last scenario, represented in Fig. [15] (c,d), it seems possible to have an axino dominated 
lightest neutralino. This is achieved with a larger value of the effective fi- term (given by Xvs) and a 
slightly lower one for the axino susy breaking parameter Mb- 

Given these results, one important issue that one would like to address concerns the modifications 
implied by our model respect to standard scenarios of neutralino densities -for instance in the MSSM 
or in the nMSSM - which require a separate investigation of the (rather large) parameter space. We 
just remark that a related analysis [30], based on an anomalous version of the MSSM which shares 
various similarities with our model, shows that for an axino-dominated LSP (light supersimmetric 
particle) - in the range between 50 GeV - 2 TeV- with a mass gap around 1-5 % between the LSP and 
the NLSP (next to lightest supersymmetric particle), the constraints from WMAP can be satisfied. 
The NLSP, in that model, has components which are typical of the (non anomalous) MSSM, with a 
dominant gaugino and/or a gaugino-higgsino projection. In the presence of extra singlets and with a 
physical axion, which is our case, this scenario should be modified even further, but we expect some 
similarities with these previous studies, especially in the neutralino sector, to hold. In a recent study 
of the axion in the MLSOM, for instance, the possibility of having the axion as a long lived particle 
require a very small mass for this particle (~ 10~^ eV) [411 . In the USSM-A the presence of an axion 
in the bosonic sector and of a neutralino in the fermionic sector as possible dark matter components 
raises the issue of the interplay between the two sectors. At the same time, in the fermionic neutral 
sector, the role of the co-annihilation becomes crucial, especially in the presence of mass degeneracy, 
which modifies substantially the neutralino relic densities already in this sector. We hope to return 
with a complete analysis of these points in the near future ^7] 
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Table 2: Parameters for the neutralino eigenvalues analysis for the charge assignment Bh^ = —3/(2^10) and 
Bh, - -l/(\/IO). 



37 



1000 2000 3000 4000 5000 

M„ (GeVl 

(a) Eigenvalues as a function of Mgt 



le-05 



500 1000 1500 2000 2500 3000 3500 4000 4500 5000 
Ms, [GeVl 



(b) Squared components of the lightest neutralino 



Figure 12: Study of the neutralino eigenvalues as a function of Stiickelberg mass Mst- 
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Figure 13: The same as Fig.([T^ but as a function of Mb and A/b- 
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Figure 14: The same as Fig. (fT2|) but as a function of My and M„ 
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Figure 15: The same as Fig. (fT2|) but as a function of gs, tan/? and wg. 
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15 Unitarity bound of the model 



Being the theory an effective description of an anomalous Lagrangean in which the presence of the axion 
is the low energy signature of a more complicated mechanism of cancellation which would eventually 
induce higher derivative terms in the effective action, it is necessary at this stage to comment about 
the unitarity of this class of models. This point has been raised in |37j and further developed in [38]. 
One of the most natural contexts for discussing unitarity is related to 2 ^ 2 processes mediated by 
BIM (Bouchiat - Iliopoulos - Meyer) amplitudes, in particular those involving gluons and photons. 
These processes exhibit an anomalous behavior when the gg — > 77 amplitude is mediated by the 
exchange in the s-channel of neutral gauge bosons that couple to the fermion loops via axial-vector 
interactions. As shown in these previous analysis, this class of amplitudes, at partonic level, violate the 
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Figure 16: BIM amplitude for gg 77 plus the amplitude obtained by the exchange of x- 



Froissart bound in the ultraviolet limit. As a matter of fact, although the Wess-Zumino counterterms 
are introduced in the Lagrangean as dimension-5 local operators to ensure the BRST invariance of the 
effective action, their contributions to the amplitudes are not sufficient to cancel the divergent behavior 
of the anomalous poles which affect the BIM amplitude shown in Fig. [16] (a). In the super symmetric 
generalization of the model that we have presented, this issue of unitarity remains basically the same 
as for the non-super symmetric case. 

As we have discussed above, in the latter case the physical axion appears as a massive degree 
of freedom in the CP-odd sector, due to the presence of a Peccei-Quinn breaking term in the scalar 
potential. After EWSB the Stiickelberg axion b is rotated directly on the physical axion x and on the 
two goldstones Gz, Gz'- Therefore, if we choose the unitary gauge, the only diagram that we can draw 
in order to erase the bad high energy behaviour of Fig. [16] (a) is the second graph (b), where the same 
amplitude of (a) is mediated by the exchange of the massive axi-Higgs, x- One can show by a direct 
study of these two graphs that there is no cancellation of these two contributions at high energy [37] . 
The problem remains also in the case of the USSM-A model discussed here. We have again a unitarity 
bound in the supersymmetric case since the only difference with respect to the non-supersymmetric 
case is the contribution of extra fermions circulating in the loops of the BIM amplitude, in particular 
the charginos. 
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16 Conclusions and Perspectives 



We have presented a generalization of the USSM in the presence of an anomalous U{1) and of a 
physical axion in the CP-odd scalar sector of the theory, model that we call the USSM-A. This model, 
which is a direct generalization of a similar construction based on the potential of the MSSM [10] . 
allows higgs-axion mixing. Both constructions are extensions of a non-super symmetric formulation, 
studied previously [S] (the MLSOM) developed in the context of orientifold vacua of string theory. 
In the case of the MLSOM, Higgs-axion mixing has been obtained by requiring that the anomalous 
gauge boson becomes massive by a combination of the Higgs and of the Stiickelberg mechanisms, with 
an axion that is part of the scalar potential. Moving to the supersymmetric case, the generalization 
of this construction - obtained by using the MSSM superpotential with an extra anomalous U{1) - is 
found to be characterized by an axino in the spectrum, which appears as a component of the neutralino 
sector, but not by an axion, since the Stiickelberg field does not acquire an axion-like coupling and 
remains a goldstone mode. The failure of the MSSM superpotential to provide such a mixing has to 
be attributed to the structure of the scalar potential of the model. Supersymmetry prohibits a term 
with a direct presence of the axion in the scalar potential, which otherwise would allow such a mixing. 

In our model the mixing occurs indirectly, but the CP-odd sector has to be non-minimal, with 
an extra singlet which is charged under the anomalous U{1). This approach, as we have emphasized, 
is quite generic, since its essential working requirement, respect to the MSSM, is the enlargement 
of the CP-odd sector with one extra SM singlet. Given these minimal requirements, which can 
be easily satisfied in rather different string vacua, these low energy effective theories capture the 
essential physical implications of several high energy scenarios, either with a low scale string scale or 
a much higher scale, as in the heterotic case. Explicit formulations of superpotentials, such as those, 
for instance, derived from free fermionic models [33], offer the natural ground where to apply the 
methodology discussed in this work. 

Anomalous C/(l)'s are quite common in string theory but can also be generated, in the correspond- 
ing effective lagrangean, by the decoupling of heavy fermions (and gauge bosons) in grand unified sce- 
narios [41j. It is then natural to ask what is left at low energy if such decoupling has indeed occurred 
at some higher scale and it reasonable to foresee that the axion is likely to play a fundamental role 
[Hj in formulating the answer to this question. Clearly, there are corrections to the action discussed 
in this work, which should be characterized by higher derivative contributions (of dimension larger 
than 5), i.e. beyond the typical Wess-Zumino terms. Arguments in favor of a possible generalization 
in this direction of the construction presented in this work have been discussed in previous works |42j 
and especially in [IT]; they are motivated by the fact that anomalies cannot be canceled with local 
counterterms. 

A related issue concerns the size of the mass of the extra Z' in the various models. It is clear 
that if its decoupling occurs at the Planck scale, then the Stiickelberg mass term takes approximately 
the value of that decoupling scale. This implies that the axion-like couplings induced at low energy 
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are also heavily suppressed. Other interactions, however, in the non-supersymmetric case, have been 
found to remain sizeable |41j . 

A final comment concerns supersymmetry breaking, which may induce phase-dependent terms in 
the potential. As discussed in [8j for the MLSOM, the axion, in that specific case, gets a sizeable 
mass which can be as large as the electroweak scale. Similar considerations could remain true in 
the supersymmetric model that we have presented, although here we have analyzed - by a deliberate 
choice - the case of a light axion, since we consider this scenario more interesting phenomenologically. 
In the presence of these phases the pseudoscalar, however, becomes massive. For instance, a mass 
region of few GeV's is certainly not excluded, as well as a scenario characterized by a very light axion 
(~ 10"^ eV), and both can be easily included within our analysis. In particular, for an axion in the 
GeV mass range, for instance, the interactions of this particle are rather similar to those of a light 
CP-odd Higgs boson, but now with extra interaction with the gauge fields, due to the anomaly, which 
are not allowed for the rest of the CP-odd sector. 
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17 Appendix A: Notations 

In this appendix we specify our notations. 
The covariant derivatives are given by 

Da = -dA - ^^""^A^/^ DA = dA + (107) 
The left /right chiral superfields in terms of field components are given in a generic form as follows 

^l{x, e, 9) = A{x) + ida^'ed^Aix) - ^e999aA{x) + V2ei^{x) 

+^99a^'9^|,^{x) + e9F{x), (108) 
v2 

^'^j^{x,9,9) = A*{x) -i9a''9d^A*{x) -^9999nA*{x) + V29^Ij{x) 

— '-^999a''di,i;{x) + 99F*{x). (109) 
v2 

A generic scalar superfield V in the Wess-Zumino gauge is given by 

V{x,9,9) = 9(T^'9[V^{x)-d^,B{x)]+999\{x)+999X{x) + 9999d{x) (110) 
where B[x) is a generic real valued scalar field. The generic expressions for the field-strengths are 
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Superfield 


Bosonic 


Fermionic 


Auxiliary 
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Q{x) 


Fq{x) 
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H.2{x) 


Mx) 


FhAx) 
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\b{x),\b{x) 


Db{x) 
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Alix) 


\y{x),\y{x) 


Dy{x) 


W^{x,9,9) 


w;(x) 


X\Yi {x), \y[ri {x) 


F>wi{x) 


G''{x,9,9) 






Dgo-{x) 



Table 3: Superfields and their components. 
852 

Qa = -^DDe'^y^^D^e^s^^ (111) 

where we have used W = t^W^ with r' being the SU{2) generators, while G = T"-G^ with T" being 
the SU (3) generators. The non supersymmetric field-strength are defined as 

pY __ aY _^ aY 

if, = d^B,-d,B^, 

G% = d^Gl-d,G",-9sr''G'Gt (112) 



Appendix B: The USSM Lagrangean 

For completeness we introduce in what follows the USSM Lagrangean that is a part of the total 
Lagrangean given by Crot = ^^USSM + l^axion + l^cs- 
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^USSM — £^lep + ^quark + ^Higgs + ^gauge + ^SMT + ^GMT 



(113) 



^quark 



iggs 



■'gauge 



l\^2g2W+gYY+gBBl ^^gzW+gyY+gBB ^ 



(114) 



= / d'd 



Q\^2gsG+2g2W+gYy+9BBQ _|_ fj\^^2g,G+gYY+gBBjj^ _|_ ^t^g2g,G+gyy+gsB^^ 



(115) 



H\e^92W+gYY+gBB fy^ ^ J^tg2g2W+ffy Y+sbS^^ + ^^e^^^-^^ + W<52(^) + Wd'^iO) 

(116) 

^ y" d^e [g^Ga + w^Wa + w^"w^ + w'^'^w^] 5'^{e) + h.c. (ii?) 

+m\HlHi + mlHlH2 + m|5^5 + (oA^i^i • if2 + /i.e.) + (oei^i • -L^ + h.c.) 
+{adHi ■ QDr + h.c.) + {auH2 ■ QUr + h.c.)] (118) 



MT 



- [MoG'^go, + M^W^Wa + MyW^'^W^ + MbW^^W^) + h.c 



5^ (9,6) 



Appendix C: The matrix 



(119) 



0^^ 



o\2 



^22 



V2VS 



V1V2 + VgV^ 
VlVS 



V\V2 + VgV"^ 
V1V2 



0, 



viifi - 2BH^gBXB + \jfl + 
V2{h + ^Bff^gBXB + ^fl + ^.g'^xl) 



2XB 



h 
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•^32 



01 



<^43 



04^4 
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BsgBvs 



2M. 



h 



viifi - 2Bh,9bXb - \Jfl + ^9 
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2Mst 



2MstVivl 
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2MstV2vl 

{vfvi + v^vl)[Blgl{vlvl + ^;2^;2) + 4M2 ^;2] ' 

2Mstvsv'^ 

BsgB\Jvlvl + v'^vl 



Blgl{vlvl + v^vl)+AMlv^ 



(120) 



[1] M. Ahlers, H. Gies, J. Jaeckel, J. Redondo and A. Ringwald , Phys. Rev. D77 (2008) 095001, 
arXiv:0711.4991[hep-ph]. 

[2] M. Ahlers, H. Gies,J. Jaeckel, J. Redondo, nd A. Ringwald , Phys. Rev. D76 (2007) 115005, 
arXiv:0706.2836[hep-ph]. 

[3] A. De Angelis, O. Mansutti and M. Roncadelli, Phys. Lett. B659 (2008) 847, 
arXiv:0707.2695[astro-ph] . 



[4] A. De Angehs, O. Mansutti, M. Persic and M. Roncadelh, (2008), arXiv:0807.4246[astro-ph]. 



46 



[5] A. De Angelis, O. Mansutti and M. Roncadelli, Phys. Rev. D76 (2007) 121301, 
arXiv:0707.4312[astro-ph] . 

[6] Z. Berezhiani, L. Gianfagna and M. Giannotti, Phys. Lett. B500 (2001) 286, hep-ph/0009290. 

[7] M.I.Vysotsky, Ya.B.Zeldovich, M.Yu.Khlopov and V.M.Chechetkin, Pisma Zh.Eksp.Teor.Fiz. 
(1978) V.27, PP. 533-536, [English translation: JETP Lett. (1978) V.27, no.9, PP. 502-505]; 

Z.G.Berezhiani, M.Yu.Khlopov and R.R.Khomeriki, Yadernaya Fizika (1990) V.52, PP. 104-109. 
[English translation: Sov.J.Nucl.Phys. (1990) V. 52, PP. 65-68]; 

Z.G.Berezhiani, A.S.Sakharov and M.Yu.Khlopov, Yadernaya Fizika (1992) V.55, PP. 1918-1933. 
[English translation: Sov.J.Nucl.Phys. (1992) V.55, PP. 1063-1071] 

[8] C. Coriano, N. Irges and E. Kiritsis, Nucl. Phys. B746 (2006) 77, hep-ph/0510332. 

[9] E. Kiritsis, Fortsch. Phys. 52 (2004) 200, hep-th/0310001. 

[10] P. Anastasopoulos, F. Fucito, A. Lionetto, G. Pradisi, A. Racioppi, Y. S. Stanev, Phys. Rev. D78 
(2008) 085014, arXiv:0804.1156[hep-th]. 

[11] J. De Rydt, J. Rosseel, Torsten T. Schmidt, A. Van Proeyen and M. Zagermann, Class. Quant. 
Grav. 24 (2007) 5201, arXiv:0705.4216[hep-th]. 

[12] M. Zagermann, (2008), arXiv:0801.1666[hep-th]. 

[13] R. Armillis, C. Coriano, M. Guzzi and S. MorelH, Nucl. Phys. B 814 (2009) 15679, 
arXiv:0809.3772[hep-ph]. 

[14] B. Kors and P. Nath, JHEP 12 (2004) 005, hep-ph/0406167. 

[15] C. Coriano, M. Guzzi, N. Irges and A. Mariano Phys.Lett.B 671 (2009) 87, arXiv:0811.0117[hep- 
ph]. 

[16] J. E. Kim and H. P. Nilles, Phys. Lett. B 138 (1984) 150. 

[17] V. Jain and R. Shrock, Phys. Lett. B 352 (1995) 83, arXiv:hep-ph/9412367. 

[18] Y. Mr, Phys. Lett. B 354 (1995) 107 [arXiv:hep-ph/9504312]. 

[19] D. Suematsu and Y. Yamagishi, Int. J. Mod. Phys. A 10 (1995) 4521 [arXiv:hep-ph/9411239]. 

[20] M. Cvetic, Durmus A. Demir, J.R. Espinosa, L.L. Everett and P. Langacker, Phys. Rev. D56 
(1997) 2861, hep-ph/9703317. 

[21] M. Cvetic and P. Langacker, Phys. Rev. D 54 (1996) 3570 [arXiv:hep-ph/9511378]. 



47 



[22] C. Balazs, M. S. Carena, A. Preitas and C. E. M. Wagner, JHEP 0706 (2007) 066 [arXiv:0705.0431 
[hep-ph]]. 

[23] U. Ellwanger and C. Hugonie, arXiv:hep-ph/0006222. 

[24] P. Anastasopoulos, M. Bianchi, E. Dudas and E. Kiritsis, JHEP 11 (2006) 057, hep-th/0605225. 

[25] D. Feldman, Z. Liu and P. Nath, JHEP 11 (2006) 007, hep-ph/0606294. 

[26] D. Feldman, Z. Liu and P. Nath, Phys. Rev. D75 (2007) 115001, hep-ph/0702123. 

[27] C. Coriano M. Guzzi and A. Mariano, in preparation. 

[28] W. Fischler, H. P. Nilles, J. Polchinski, S. Raby and L. Susskind, Phys. Rev. Lett. 47 (1981) 757 

[29] J. Preskill, Ann. Phys. 210 (1991) 323. 

[30] C. Coriano, N. Irges and S. MorelH, JHEP 07 (2007) 008, hep-ph/0701010. 

[31] C. Panagiotakopoulos and K. Tamvakis, Phys. Lett. B446 (1999) 224, hep-ph/9809475. 

[32] S.A. Abel, S. Sarkar and P.L. White, Nucl. Phys. B454 (1995) 663, hep-ph/9506359. 

[33] D.J. Miller, R. Nevzorov and P.M. Zerwas, Nucl. Phys. B681 (2004) 3, hep-ph/0304049. 

[34] C. Coriano, N. Irges and S. Morelh, Nucl. Phys. B789 (2008) 133, hep-ph/0703127. 

[35] C. Amsler et al., Physics Letters B667, 1 (2008). 

[36] R. Armillis, C. Coriano and M. Guzzi, JHEP 05 (2008) 015, arXiv:0711.3424[hep-ph]. 

[37] C. Coriano, M. Guzzi and S. Morelh, Eur. Phys. J. C55 (2008) 629, arXiv:0801.2949[hep-ph]. 

[38] R. Armillis, C. Coriano, M. Guzzi and S. Morelh JHEP 0810 (2008) 034, arXiv:0808.1882[hep-ph]. 

[39] J. Kalinowski, S.F. King and J.P. Roberts, JHEP 0901 (2009) 066, arXiv:0811.2204[hep-ph]. 

[40] A. Lionetto and A. Racioppi, arXiv:0905. 4607 [hep-ph]; F. Fucito, A. Lionetto, A. Mammarella 
and A. Racioppi, arXiv:0811.1953[hep-ph]. 

[41] C. Coriano and M. Guzzi, arXiv:0905.4462[hep-ph] . 

[42] C. Coriano and N. Irges, Phys.Lett.B651 (2007), 298, hep-ph/0612140. 

[43] A.E. Faraggi, E. Manno and C. Timirgaziu, Eur. Phys. J. C50 (2007) 701, hep-th/0610118. 



48 



